RM E57K19a 


RESEARCH MEMORANDUM 

CORRELATION OF TURBULENT HEAT TRANSFER IN A TUBE FOR 

THE DISSOCIATING SYSTEM N 2 0 4 -=* 2N0 2 

* 

By Richard S. Brokaw 

Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

WASHINGTON 

March 5, 1958 


NACA EM E57KL9a 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
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CORRELATION OF TURBULENT HEAT TRANSFER IN A TUBE 
FOR THE DISSOCIATING SYSTEM N 2 0 4 * 2N0 2 
By Richard S. Brokaw 

SUMMARY 

Experimental data on turbulent heat transfer in a tube for the dis- 
sociating system N 2 0 4 2N0 2 are in accord with the usual Nusselt- 

Prandtl-Reynolds number correlation for this type of convective heat 
transfer. 

The thermal conductivities and kinematic viscosities needed to 
formulate the Nusselt, Reynolds, and Prandtl numbers have been computed 
using rigorous expressions from the kinetic theory of gases. The thermal 
conductivity includes a large contribution arising from the diffusional 
transport of chemical enthalpy. 

Three methods of selecting average property values for the boundary 
layer were tested. All gave good correlation using the equilibrium trans- 
port properties. 

An approximate method using frozen gas properties gave somewhat poor 
results, but should be useful for complex systems where exact computation, 
would be prohibitively tedious. This approximate method should also be 
applicable to systems in which chemical equilibrium is not attained, pro- 
vided the enthalpy of the gases in the free stream and at the wall is 
known. 


INTRODUCTION 

Currently there is considerable interest in heat transfer in disso- 
ciating gases, because such gases may exist in the exhaust nozzles of jet 
engines or within the boundary layer on hypersonic aircraft. In reacting 
gases, heat transport may be considerably higher than in "frozen" (non- 
reacting) mixtures. Large amounts of heat may be carried as chemical 
enthalpy of molecules which diffuse because of concentration gradients. 
These gradients arise because the gas composition varies with temperature. 
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For example, in a gas which absorbs heat by dissociating as the temper- 
ature is raised, heat is transferred when a molecule dissociates in the 
high- temperature region and diffuses toward the low- temperature region 
(there is a lower concentration of dissociated molecules at low temper- 
ature) . In the low- temperature region the gas reassociates, releasing 
the heat absorbed from the high- temperature dissociation. 

Heat transfer in practical systems with chemical reaction has been 
theoretically examined (e.g., refs. 1 and 2). Parallel theoretical 
studies of the thermal conductivity of reacting gases have been made : 
in reference 3 a general expression for the thermal conductivity of gas 
mixtures in chemical equilibrium is developed, while in reference 4 
systems with finite reaction rates are considered. 

Recently Beal and Lyerly (ref. 5) have experimentally investigated 
turbulent heat transfer in a tube for the dissociating system N 2 O 4 2N0 

Essential features of their apparatus included a 3/4-inch-diameter test 
section with a rounded entrance. Heat was added at a thermally guarded 
hot spot 23 diameters from the tube entry. The guard heater was short, 
so that the length- diameter ratio for the thermal boundary layer was only 
2. Tests were carried out with air in the Reynolds number range from 
2500 to 34,000 and with N 0 O 4 2N0-? mixtures in the range from 18,000 
to 104,000. 

Beal and Lyerly were unable to analyze their results fully, because 
of the dearth of experimental conductivity and viscosity data for the 
w 2°4 * 2N0 2 s y stem in the literature. This report shows that the data 

of Beal and Lyerly are in complete accord with the usual Nusselt-Prandtl- 
Reynolds number correlations for this type of convective heat transfer. 
Thermal conductivities and viscosities for the N 2 O 4 ^ 2 NO 2 system 

(needed to formulate the Nusselt, Prandtl, and Reynolds numbers) have 
computed using rigorous kinetic theory expressions. 


SYMBOLS 

The following symbols are used in this report: 
dimensionless constant, eq. (l) 
specific heat at constant pressure 
binary diffusion coefficient 
tube diameter 
mass velocity 


been 


A 

D 

d 

G 
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Ah specific enthalpy difference between wall and fluid bulk 

h heat- transfer coefficient 

k thermal conductivity 

Nu Kusselt number, hd/k 

Pr Prandtl number, c^p/k 

Re Reynolds number, Gd/p or Gd/pv 

T absolute temperature 

8 pDCp,f/kf (dimensionless) 

p absolute viscosity 

v kinematic viscosity 

p density 

<p generalized temperature dependent property 

Subscripts : 

B bulk fluid 

e chemical equilibrium 

f chemically frozen 

i contribution due to internal energy 

m value for mixture of monatomic gases 

r contribution due to chemical reaction 

w wall 

Superscripts : 

— mean value 

* evaluated at reference condition 
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METHODS OF ANALYSIS 
Transport Properties 

The computed thermal conductivities, kinematic viscosities, and 
Prandtl numbers for the NgO^ 2N0 2 system at 1 atmosphere are shown 

in figure 1. The methods of computation are detailed in the appendix. 
Calculated conductivities agree with recent experimental determinations 
(ref. 6), which are also plotted in figure 1. Frozen thermal conduc- 
tivities and Prandtl numbers are shown by the dashed lines. (The frozen 
values refer to equilibrium gas composition without chemical reaction.) 

The horizontal bars across the top of the figure indicate the temperature 
differences between the wall and bulk fluid for the individual N 2 04 ^ 2 N 0 2 

heat-transfer tests of Beal and Lyerly; approximate Reynolds numbers are 
given adjacent to the bars. (The exact values of the Reynolds numbers 
depend on the way mean property values are chosen.) Over the tempera- 
ture range of the heat-transfer experiments, the conductivity varies by 
a factor of 2.2, the kinematic viscosity by a factor of 2.3, and the 
Prandtl number by a factor of 1.4. These variations are much larger than 
those encountered with nonreacting gases; over a comparable temperature 
range, conductivity might increase by 10 to 15 percent and kinematic vis- 
cosity by 30 to 40 percent, while Prandtl number would remain substan- 
tially constant. Accordingly, for the NgO^ 2N0 2 system, suitable 
mean property values for correlation must be selected. 


Selection of Mean Property Values 

In order to formulate the Reynolds, Nusselt, and Prandtl numbers, 
mean values of viscosity, thermal conductivity, and Prandtl number are 
required. In actual fact, a mean value of kinematic viscosity is_re- 
quired in the Reynolds number, which may be written Gd/p = Gd/p-g v ; 
hence, 


Mean properties obtained by the following three methods were tested: 


(l) Temperature- integrated properties. If a property is designated 
as <*, then the integrated average value is 
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This type of average suggests itself naturally for thermal conductivity 
because it is the proper average for heat transfer in the absence of con- 
vection (e.g., thermal- conductivity cells). 

(2) Reference- temperature properties: Properties are taken at a 

reference temperature , as suggested in reference 7 . 

T* = T b + 0.58 (T w - T b ) 

(3) Reference- enthalpy properties: Properties are evaluated at a 

temperature corresponding to a reference enthalpy h*, defined as 

h* * h B + 0.58 (l^ - h B ) 


COMPARISON OF CORRELATIONS 

Correlations Using Equilibrium Transport Properties 

All three methods of averaging properties give good correlation of 
the air and N 2 0 4 £ 2N0 2 heat- transfer data of Beal and Lyerly. Figure 

2 shows a log-log plot of Nu/Pr 1 / 3 against Reynolds number for 
temperature-integrated average properties. A mean kinematic viscosity 
is required in formulating the Reynolds number. If a mean value of the 
absolute viscosity is used, correlation is very poor, with errors of up 
to fifty percent in the Reynolds number. This is shown in figure 3, a 
reference- enthalpy equilibrium- property correlation, in which a mean 
absolute viscosity was used. The experimental data for each of the three 
correlations using mean kinematic viscosity were fitted by an equation 
of the form 

Nu/Pr 1 / 3 = ARe 0,79 (l) 

or 

log Nu - - log Pr = log A + 0.79 log Re (la) 

3 

For each correlation, the constant A was chosen to give a least- squares 
fit of all the data (air and N 2 0 4 ) at Reynolds numbers greater than 

9000 to equation (la) . (The air points at Re < 9000 were excluded be- 
cause they lie in the transition region.) Results of the correlations 
are summarized in the first three columns of table I 1 . 


lit might be noted that in the Reynolds number range under consider- 
ation (9000 to 110,000) equation (l), with the A values of table I, 
agrees to within 4 percent with Colburn's correlation (ref. 8): 

Nu/Pr 1 / 3 = 0.023 Re 0,8 . This is somewhat surprising, in view of the 
short length-diameter ratio for the thermal boundary layer in Beal and 
Lyerly' s apparatus. 
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From table I it is seen that all three correlations represent the 
N 2 0 4 * 2N0 2 data very well. (The standard deviations for the air data 

and over-all correlations merely indicate the consistency of the corre- 
lations.) For engineering heat- transfer calculations the reference- 
enthalpy or reference- temperature methods would ordinarily be used, be- 
cause with these methods, transport properties need be calculated at only 
one condition. 

However, table I shows that the integrated property correlation is 
slightly better. This is merely a tentative conclusion,' however, since 
the differences between the correlations are largely the result of a 
single data point (albeit the point with the largest temperature differ- 
ence). More experimental measurements with larger temperature differ- 
ences are required to firmly establish the most reliable method for com- 
puting mean properties. In any event, for all three correlation methods, 
deviations for N 2 0 4 and air data points are similar. This fact strongly 
supports the notion that heat-transfer correlations for nonreacting sys- 
tems may be applied to dissociating systems, provided that transport 
properties are properly computed and suitable mean values are chosen. 


Approximate Correlation Using Frozen Transport Properties 

In reference 3 it is shown that the equilibrium and frozen thermal 
conductivities may be related to the equilibrium and frozen specific heats 
by the equation 



If the dimensionless group 6 is assumed to be unity, then 

k e s k f c p,e/ c p,f • (Aerodynamicists sometimes refer to 8 as the "Lewis 

number." However, this usage is not generally accepted.) Hence, for a 
system with a finite temperature difference, it seems reasonable to de- 
fine an approximate mean thermal conductivity by 



with the frozen specific heat and conductivity evaluated at a reference 
temperature. Subject to this approximation, the Nusselt and Prandtl num- 
bers are 

c. fAT 
NUe B Nu f -te— 


Pr e £ Pr f 
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A reference- temperature correlation using frozen properties is shown 
in figure 4. (if the correction term Cp^jAT/Ah were not applied to the 

frozen Nusselt number, Nusselt numbers for the ^ 2W0g system would 

be eight- to ninefold too high.) As expected, the correlation of figure 
4 is poorer than the correlation of figure 2; this is reflected in a 
higher standard deviation, shown in the fourth column of table I. In 
reality the deviations are systematic with 6 and might amount to 10 
percent or more under some conditions. For this system 6 actually 
ranges from 1.13 at 310° K to 0.64 at 375° K. The frozen property cor- 
relation can be improved by introducing a correction term involving 5 
and plotting 


(NUf/Prl/^Ccp^/Ah)^! + (6 - 1)^1 - 


-2/3 


against Reynolds number. The correction term is reminiscent of similar 
corrections used in treatments of hypersonic stagnation- point heat trans- 
fer (refs. 2 and 9). However, heat transfer in multicomponent mixtures 
cannot, in general, be characterized by a single diffusion coefficient, 
and thus the correction term is not generally useful^. 


The approximate correlation without the correction should be appli- 
cable to many complex systems where an exact computation would be pro- 
hibitively tedious (e.g., jet-engine and rocket exhausts). Indeed, it 
should be applicable to systems in which chemical equilibrium is not at- 
tained, provided the chemical compositions (and hence, enthalpy) of the 
gases in the free stream and at the wall are known. 


CONCLUSIONS 

Experimental data for turbulent heat transfer in a tube have been 
analyzed and correlated. The following conclusions may be drawn: 

1. Experimental heat-transfer data in a tube for the dissociating 
system N2O4 ^ 2NO2 correlate well with air data using transport prop- 
erties computed on the basis of chemical equilibrium. 

2. In computing Reynolds number a mean kinematic viscosity is 
required. 


^For multicomponent mixtures & may be defined as 
[(k e /kf) - 1 ]/[ ( c p,e/ c p,f ) “ and the correction term is valid. This 

is not useful, however, since the equilibrium properties calculated to 
obtain & might better be applied directly in the correlations of the 
preceding section. 
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3. The best correlation was obtained using temperature- integrated 
properties. Reference- temperature and reference- enthalpy properties were 
almost as good. 

4. Errors are somewhat larger for an approximate correlation which 
uses frozen properties^ In this method the frozen Nusselt number is cor- 
rected by the factor Cp^^AT/Ah and it is assumed that pDCp^p/kp = 1. 

(Mean value of the specific heat at constant pressure and chemically 
frozen conditions, Cp^p; temperature difference between wall and fluid 

bulk, AT; specific enthalpy difference between wall and fluid bulk, Ah; 
density, p; binary diffusion coefficient, D; specific heat at constant 
pressure and chemically frozen conditions, Cp^p; chemically frozen ther- 
mal conductivity, kp.) 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, December 6, 1957 
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APPENDIX - CALCULATION OF EQUILIBRIUM PROPERTIES 
FOR THE N 2 0 4 ^ 2N0 2 SYSTEM 

Transport-property calculations were carried out assuming a Lennard- 
Jones (6-12) potential for the molecular interactions (ref. 10) . Gas 
compositions were obtained from the equilibrium data of Bodenstein (ref. 
ll). The means by which various properties were computed are as follows: 

Force constants . - In order to compute the transport properties, 
force constants a and e/k characteristic of molecular interactions 
must be estimated. For the N 2 0 4 molecule, the force constants were es- 
timated from critical temperature and density (ref. 10, eqs. (4.1-17) 
and (4.1-18)). The critical temperature and density for pure N 2 0 4 were 

in turn estimated from the boiling point and density of liquid N 2 0 4 by 

the methods of references 12 and 13. (it was assumed that the boiling 
point and liquid density of N 2 0 4 2N0 2 equilibrium mixtures is close 

to the value which would be expected for pure N 2 0 4 .) Force constants 
chosen for N 2 0 4 were: o = 4.74 A and s/k = 383° K. 

The force constants for the NO? molecule were assumed to be similar 

c o 

to those of C0 2 and N 2 0 (ref. 10, p. 1111 ) ; values used were a = 3.9 A 

and e/k *> 230° K. 

For the N 2 0 4 -N0 2 interaction, a was taken as the arithmetic mean 
of N0 2 and N 2 0 4 values; e/k was taken as the geometric mean (ref. 10, 
eqs. (8.4-8) and (8.4-9)). 

These force constants differ somewhat from those used in reference 
3; they are believed to be more realistic estimates. 

Viscosities . - Viscosities were calculated by equation (8.2-22) of 
reference 10. 

Thermal conductivities . - The thermal conductivity of a reacting 
gas mixture may be broken into three portions: 


k e - ^m + *1 + *r ( 3 ) 

The first term on the right side of equation (3) represents collisional 
kineti c . energy transport; it was computed using equation (8.2-36) of ref- 
erence 10. The second term (in the nature of an Eucken correction) rep- 
resents a diffusional transport of internal energy (rotation and vibra- 
tion) and was calculated using equation (75) of reference 14.. The third 
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and dominant term for these dissociating mixtures is the increase in con- 
ductivity due to chemical reaction. It was computed from equation (l4) 
of reference 3. The frozen thermal conductivity is simply the first two 
terms of equation (3 ) . 

Diffusion coefficients . - Diffusion coefficients were computed 
according to equation (8.2-44) of reference 10. 

Specific heats . - The specific heat may he broken up into a frozen 
part and a contribution due to chemical reaction. The frozen specific 
heats were computed from the known vibrational frequencies of the NOg 
and N2O4 molecules. The contribution due to chemical reaction was com- 
puted according to reference 3. 

Heat of reaction . - The heat of reaction at 298.16° K was computed 
from the heats of formation of the N0 2 and N 2 0 4 molecules presented in 

reference 15. Values at other temperatures were obtained by correcting 
the 298.16° K value for the integrated heat capacity difference between 
products and reactant. 

Calculated values of viscosity, density divided by pressure, frozen 
and equilibrium thermal conductivity, and frozen and equilibrium specific 
heat are contained in table II. Values were computed at 10° intervals 
from 290° to 490° K and for pressures of 1.0, 0.74, and 0.33 atmospheres. 
(Additional thermal conductivities computed for pressures of 0.5, 0.2, 

0.1, 0.05, and 0.02 atmospheres may be found in reference 6.) Only the 
1.0- atmosphere values were used in the correlations of «. this report. 
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TABLE II. - COMPUTED PROPERTIES FOR N 2 0 4 - 2N0 2 EQUILIBRIUM 
(a) Pressure, 1 atmosphere. 


Temperature , 
°K 

Viscosity, 

poises 

Equilibrium thermal 
conductivity, 
cal/(cm)(sec) (°K) 

Frozen thermal 
conductivity. 
cal/(cm) (sec)( 6 K) 

Density + 
pressure, 
gm/(cm 3 )(atm) 

Equilibrium 
specific heat, 
cal/gm 

Frozen 
specific 
heat , 
cal/gm 

290 

1.124X10" 4 

20.083X10“ 5 

3.110X10" 5 

3.414*10' 3 

1.0168 

0.2017 

300 

1.192 

25.668 

3.337 

3.130 

1.2973 

.2042 

310 

1.270 

31.176 

3.587 

2.832 

1.6051 

.2060 

320 

1.356 

35.233 

3.860 

2.540 

1.8730 

.2075 

330 

1.445 

36.366 

4.143 

2.272 

2.0067 

.2083 

340 

1.534 

33.977 

4.427 

2.043 

1.9395 

.2089 

350 

1.615 

29.086 

4.693 

1.864 

1.6893 

.2094 

360 

1.689 

23.300 

4.938 

1.726 

1.3483 

.2099 

370 

1.754 

18.037 

5.163 

1.621 

1.0150 

.2106 

380 

1.811 

14.087 

5.366 

1.542 

.7556 

.2115 

390 

1.865 

11.262 

5.561 

1.479 

.5650 

.2125 

400 

1.915 

9.478 

5.745 1 

1.428 

.4415 • 

.2138 

410 

1.961 

8.375 

5.923 

1.385 

.3622 

.2151 

420 

2.005 

7.742 

6.098 

1.346 

.3134 

.2166 

430 

2.048 

7.397 

6.272 

1.312 

.2831 

.2181 

440 

2.091 

7.219 

6.446 

1.280 

.2634 

.2196 

450 

2.133 

7.143 

6.617 

1.250 

.2503 

.2211 

460 

2.174 

7.163 

6.791 

1.221 

.2429 

.2226 

470 

2.215 

7.223 

6.958 

1.195 

.2382 

.2240 

480 

2.254 

7.315 

7.127 

1.169 

.2353 

.2255 

490 

2.294 

7.436 

7.300 

1.145 

.2341 

.2270 


(b) Pressure, 0.74 atmosphere. 


290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 


1.134xl0" 4 

1.206 

1.288 

1.377 

1.468 

1.554 

1.633 

1.702 

1.763 

1.817 

1.870 

1.918 

1.963 

2.006 

2.049 

2.092 

2.133 

2.175 

2.215 

2.254 

2.294 


22. 332x10 " 5 
28.201 
3o.**98 
36.578 
36.073 
32.025 
26.197 
20.304 
15.542 
12.201 
9.939 
8.576 
7.769 
7.331 
7.114 
7.024 
7.008 
7.069 
7.154 
7.267 
7.401 


3.136K10" 5 

3.373 

3.633 

3.916 

4.203 

4.482 

4.740 

4.975 

5.188 

5.384 

5.573 

5.753 

5.929 

6.102 

6.275 

6.448 

6.618 

6.792 

6.959 

7.128 

7.300 


3.352X10"’ 

3.053 

2. 745 

2.449 

2.189 

1.974 

1.811 

1.688 

1.596 

1.526 

1.469 

1.422 

1.380 

1.343 

1.310 

1.278 

1.249 

1.221 

1.194 

1.169 

1.145 


1.1407 

0.2015 

1.4500 

.2039 

1.7685 

.2056 

2 . 0051 

.206 9 

2.0543 

.2076 

1.8776 

.2081 

1.5482 

.2086 

1.1803 

.2092 

.8677 

.2101 

.6410 

.2111 

.4841 

.2123 

.3866 

.2136 

.3258 

.2150 

.2890 

.2164 

.2666 

.2180 

.2522 

.2196 

.2427 

.2210 

.2377 

.2226 

.2344 

.2240 

.2328 

.2255 

.2322 

.2270 


(c) Pressure, 0.33 atmosphere. 


290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 


1.166x10 

1.250 

1.342 

1.435 

1.523 

1.600 

1.667 

1.727 

1.780 

1.828 

1.877 

1.922 

1.966 

2.009 

2.051 

2.093 

2.134 

2.175 

2.215 

2.255 

2.294 


29.106x10 

34.829 

37.844 

36.481 

31.189 

24.217 

18.109 

13.554 

10.557 

8.752 

7.663 

7.078 

6.781 

6.667 

6.657 

6.711 

6.796 

6.918 

7.048 

7.191 

7.346 


3.221x10 

3.488 

3.776 

4.072 

4.352 

4.608 

4.835 

5.042 

5.234 

5.414 

5.593 

5.766 

5.937 

6.108 

6.279 

6.451 

6.620 

6.794 

6.960 

7.129 

7.301 


3.152x10 

2.819 

2.497 

2.216 

1.994 

1.828 

1.709 

1.621 

1.554 

1.499 

1.452 

1.411 

1.373 

1.339 

1.307 

1.276 

1.247 

1.220 

1.194 

1.169 

1.145 


1.5476 

0.2010 

1.9064 

.2030 

2.1643 

.2043 

2.1774 

.2051 

1.9129 

.2057 

1.4900 

.2063 

1.0850 

.2070 

.7694 

.2080 

.5566 

.2092 

.4248 

.2105 

.3422 

.2118 

.2943 

.2133 

.2650 

.2148 

.2494 

.2163 

.2398 

.2179 

.2342 

.2195 

.2307 

.2210 

.2293 

.2225 

.2287 

.2239 

.2287 

.2254 

.2293 

.2270 






Prandtl number. Kinematic viscosity. Thermal conductivit; 

Pr - v, cm 2 /sec cal/(cm) (sec ) (°1 
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Figure 2. - Correlation of data with temperature- integrated equilibrium properties. 
(Mean kinematic viscosity used to formulate Reynolds number.) 
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Figure 3. - Correlation of data with reference-enthalpy equilibrium properties. 
(Mean absolute viscosity used to formulate Reynolds number. ) 
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